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Abstract 

We introduce a new technique to calculate perturbative corrections to neutron-deuteron (nd) 
scattering that does not require calculation of the full off-shell scattering amplitude. Its relation 
to the more familiar partial-resummation technique is explained. Also included is a calculation of 
the SD-mixing term that occurs at next-to-next-to-leading-order (NNLO) in pionless effective field 
theory (EFTy). Using the new technique with the SD-mixing term a complete strictly perturbative 
phase-shift analysis of nd scattering is performed up to NNLO including eigenphases and mixing 
angles. This is compared to potential model calculations and good agreement is found with the 
eigenphases and some of the mixing angles at low energies. 
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I. Introduction 



At sufficiently low energies (E < m^/Mjv), where pions can be integrated out, nuclear 
physics can be described by pionless effective field theory(EFT ? /) in which only nucleon de- 
grees of freedom as well as possible external currents appear. (For a review of EFTy see jl].) 
This approach has been used quite successfully to calculate nucleon-nucleon (NN) scattering 
2-p], electromagnetic form factors of the deuteron, and the neutron proton capture process 
6, 7] in the two-body sector. It has also been used to study iViV parity- violation js, 9] and 
neutrino-deuteron processes [loj]. In the three-body sector calculations have been carried out 
for the bound state properties of Helium-3 parity- violation in neutron-deuteron (nd) 



interactions 



12 



19 



131 ] . and nd scattering 14j-ll8| as well as proton deuteron (pd) scattering 



20| in which Coulomb forces were treated perturbatively. However, in the three-body 
case for nd scattering no strictly perturbative calculation has been carried out to NNLO 
since it apparently requires the full off shell leading-order (LO) scattering amplitude, which 
is numerically expensive. In order to calculate nd and pd scattering to higher orders the 



partially resummed approach 16[ has been used in which certain classes of diagrams are 
summed to all orders. The calculation contains all the correct diagrams to the order one 
is working, but also contains certain higher order diagrams. This method leads to good 
results for nd scattering in both the effective range expansion (ERE) parametrization and 



Z parametrization 



16 



211 ] . However, certain problems occur in the quartet S-wave channel 



above the deuteron breakup threshold (DBT), namely the imaginary part of the phase shift 
at next-to-leading-order (NLO) and NNLO is negative and hence unphysical. A possible way 
around this is to resum to all orders the effective range in the deuteron propagator and then 
use this deuteron propagator to calculate the NNLO quartet S-wave scattering amplitude 
This approach gives physical results for the phase shifts but introduces spurious poles 



in the deuteron propagator. These poles occur for momenta greater than t 



However, they will still introduce issues with the numerical solution 



re cutoff of EFT^. 



181 ] . A perturbative 



calculation using short range effective fie 
has been carried out by Ji and Phillips 



d theory (SREFT) for a trimer of Helium-4 atoms 



221 ] . They calculate the full off shell leading-order 



scattering amplitude and use it to numerically solve all necessary Feynman diagrams up to 
and including NNLO. This calculation is similar to nd scattering in the doublet S-wave chan- 
nel. In this paper we propose a new method to calculate to arbitrary order perturbatively 



without the need to calculate the full off-shell LO scattering amplitude. The method puts 
contributions into the inhomogeneous part of the integral equation that depend on ampli- 
tudes that are lower order than the order at which one is working. At each order the integral 
equation has the same kernel but a different inhomogenous term. The paper is organized 
as follows. In section II we introduce the two-body Lagrangian and then in section III we 
calculate nd scattering in the quartet channel to NNLO and introduce the new technique. 
Section IV extends this to the doublet channel and in section V the SD-mixing terms for 
nd scattering occurring at NNLO are calculated. In section VI we describe our phase shift 
analysis conventions. In section VII we analyze the results and make conclusions in section 
VIII. 

II. Two-Body Lagrangian 

The two-body Lagrangian in the auxiliary field formalism is given by 



t]N T PiN + h.c. 



(1) 



V 
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where tj (s a ) is the spin-triplet (spin- singlet) dibaryon field. The auxiliary field Lagrangian 
can be shown to be equivalent to a Lagrangian containing just nucleon fields by integrating 
out the auxiliary fields and performing a field redefinition [23[. The projector Pj = ^go^cr^ 
(P a = -^T 2 T a (j2) projects out the spin-triplet iso-singlet (spin-singlet iso-triplet) combination 
of nuclei. The term is sub-leading compared to . Thus the LO bare deuteron 

propagator is simply given by i/hSJU and at LO is dressed by an infinite number of nucleon 
bubbles as in Fig[D This bubble sum can by solved explicitly by means of a geometric series 
yielding 



where the [i dependence is obtained by using dimensional regularization with the power 
divergence subtraction scheme (PDS) [24j in which poles occurring in three dimensions are 



added back to our expressions. The coefficients are then chosen such that the deuteron pole 
is in the right location, giving 



A (-i) _ M N 



yf 47r 



(7t~/4 (3) 



where jt = 45.7025 MeV is the deuteron binding momentum. At NLO the deuteron prop- 
agator receives a correction from the deuteron kinetic term and the NLO correction A^ , 
and at NNLO it receives two such corrections as shown in FigfJJ Thus the full deuteron 
propagator up to and including NNLO is given by 

iOf MO ( P „,p) = iA to (po,p) (1 + A 10 (po,p)(a;§ +p„ - (4) 

+ (I^, W )»(A* +J *-i)' 

(3g \ (3g \ 

There exists two different approaches by which to fit the parameters A)^ A| ^ , and 
y t . In the effective range expansion (ERE) parametrization one insures that the deuteron 
pole is given correctly and that the deuteron pole residue Z t = 1/(1 — r ytPt) is given correctly 
in an expansion of the effective range about the deuteron pole, p t = 1.65 fm. The other 



approach, termed the Z-parametrization [2l|, also fits to the deuteron pole, but fits to the 
deuteron pole residue exactly at NLO. The deuteron propagator in both parametrizations 
has been used extensively throughout the literature and we merely quote the results for the 
deuteron propagator in the Z-parametrization [18], which is given by 



Ani 1 

It ~ \/? - m nPo - it 



I + ^— ^ ( lt + M NPo - ie] + (%-^\ (%--M N po-it) + 



2 lt \ V 4 \ 2 lt J \A 



NNLO 



LO NLO 

and the resulting constraints on the coefficients which are 

JL _ M n Zt-1 (3 5l) _ 2y 2 t fj-jt ( 3 5l) _ jl_ 



yj M lt 1 + (Z t - 1) ' M N Z t -V (°) ~ M 1 



N 
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(LO) 

X X 



(NLO) (NNLO) 



( 3 . 

nucleon bubbles to give the LO dressed deuteron propagator. At NLO the dressed 



FIG. 1: At LO the bare deuteron propagator i/hSJU is dressed by an infinite number of 



deuteron propagator receives one effective range correction which comes from the deuteron 

( 3 Si) 

kinetic term and the NLO correction ALs . Finally at NNLO the dressed deuteron 
propagator receives two such effective range corrections 

The form of the So dibaryon propagator is obtained analogously to that of the deuteron 
propagator. However, in the ERE parametrization one expands about zero momentum, since 
their is no physical bound state pole. Then one fits to obtain the correct scattering length. 
At higher orders in the ERE one fits to the effective range. For the Z-parametrization one 
insures that the virtual-bound state pole at LO in the channel is reproduced correctly. 
At higher orders the correct residue about the virtual-bound state pole is reproduced and 
the pole is not changed. Since there is no physical data for the virtual-bound state pole one 
must infer its location from scattering data, and doing so one finds the virtual-bound state 
binding momentum 7 S = —7.8902 MeV, and the residue about the pole Z s = .9016 
The dibaryon propagator in the Z-parametrization is given by 



18|- 



iD^° M ) = 1 (7) 

NVs Is ~ J f - M NPo - ie 



X 



z s -i ( [$ 2 ~ T\ , fz s -i\ 2 fp 2 w 



2 7 s \ ' V 4 — / V 2 7 , J V 4 



Is ~\~ \ . M NPo - ie + ^ - M NPo - Y s + ■ 



LO NLO 

and the resulting constraints on the coefficients in 



NNLO 



J_ M 2 N Z s -1 ( i So ) _ 2y 2 fi- ls ( i 5o) _t|_ 

vl M ls i + ( Zs -iy^ m n z s -v M N w 

For the sake of later convenience we adopt the notation (p , p) where n=0,l,2,.. refers 
to LO, NLO, NNLO, and higher orders respectively. Thus D[°\p ,p) is the LO deuteron 



5 



propagator and D^{p ,p) only picks out the NNLO piece of the deuteron propagator, the 
part in the brackets of Eq. fl5]) that are under-braced with NNLO times the piece outside of 
the brackets. An analogous notation is also used for the dibaryon propagator. 

Finally we note that the deuteron wavefunction renormalization is given by the residue 
of the deuteron propagator about the deuteron pole 



J D 




+ 

N"LO 



(9) 



In the Z-parametrization the deuteron pole residue is given exactly at NLO and there- 
fore there are no NNLO or higher order corrections to the deuteron wavefunction renor- 
malization. For later convenience we will denote Zlo, Z^lo, and Z^nlo as the part of 
the deuteron wavefunction renormalization occurring at each order in EFTy. In particular 

Z L o = (87T7t)/(Mjv-yf) , Z NLO = Z LO Z t , and Z NNLO = 0. 

III. Quartet channel 

In order to calculate nd scattering in the quartet channel at LO in EFTy one must sum 
an infinite number of diagrams. However, unlike the two-body case the diagrams do not 
factorize and cannot be solved analytically by means of a geometric series, rather the infinite 
sum of diagrams can be solved numerically by means of an integral equation represented in 
Fig. [2] and originally given by Skornyakov and Ter-Martirosian 25] . The integral equation 
for the half off-shell scattering amplitude represented in Fig [2] is given by Eq. fflQj) 161 ] . 






FIG. 2: Integral equation for quartet channel at LO 
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(W i ) pb (k,p,/i) = V l^a'fj h a —, J—? = + (10) 

V / aa ^ 2 V ,a a k 2 7 t 2 | h (k+p) 2 , • V 7 

4Af]v Afjv ' 2Mat 

+ f (*V)?*2 / -^W k )Z& q, & + go) x 

x iD?) f — - 4" + ft + g , cA _ * x 

* \4M N M N y °' 4 l 



2M N ^ u 2M 



9o - o^7 + « 



x 



4Mjv M N 1 1 ^ 2M N 

where i (j) is the initial (final) deuteron polarization, a (/3) the initial (final) nucleon spin, 
and a (b) the initial (final) nucleon isospin. The incoming momentum in the center of mass 
(cm.) frame is k and the outgoing momentum in the cm. frame is p. Finally the parameter 
h = (k 2 — p 2 )/2Mat measures the off-shellness of the outgoing particles. Note that when 
|p| = |k| the amplitude is full on shell and h = 0. 

We now simplify Eq. (fTOj) by integrating over the energy and picking up a simple pole 
in the complex energy plane. Then we project the spin in the quartet channel, set h = 
(k 2 — p 2 )/2Mtv, and perform a partial wave decomposition of the amplitude yielding [if]]. 



pk \ pk 

2 [ k 1 2I/I , 1 1 fp 2 + q 2 -M N E-ie 
+ - / dqq 2 t l q (k,q) =—Q l P - * ^ 



Q l (a)= 1 - fdxZQ, (12) 



where 



2 J_ 1 x + a ' 

is equal to the Legendre polynomials of the second kind up to a factor of (— 1)'. The 
energy appearing in Eq. ( ITT]) is defined as E = — (Note the subscript q on the 
amplitude t refers to the quartet channel.) The partial wave decomposition changes the 
integral equation from three dimensions to one, but there is now an infinite set of integral 
equations. However, at low energies only a few partial waves are needed to obtain sufficient 
convergence and for our purposes we calculate up to and including G-waves. To solve the 



integral equations numerically we employ the Hetherington-Schick method 



2fH2g|], in which 



one solves the integral equation along a contour rotated into the complex plane thereby 
avoiding the logarithmic singularities and the simple pole from the deuteron propagator. 
This is valid as long as no singularities exist within the contour and this is satisfied for 
our equation. Finally one can use the integral equation and the solution of the amplitude 
along the contour to solve for the amplitude along the real axis again provided there are no 



singularities within the resulting contour, which can also be shown to be true 29]. 



We note that in Eq. ( TTTj) the integral ranges from to A instead of to infinity. The 
use of a cutoff for the integral serves two purposes. First, it gives a cutoff regularization for 
potential divergences. Second, since we will have to solve this integral equation numerically 
we have to impose some cutoff anyway since we cannot integrate to infinity. One may 
be concerned that we are using two different regularization schemes here. We are using 
dimensional regularization for the deuteron propagator and a cutoff regularization for the 
integral equation. However, any effects from choosing different regularization schemes are 
of higher order 15 1 . 



At NLO the scattering amplitude is given by the diagram in Fig. [3] where a single 
insertion of the effective range correction is put between two half off-shell LO scattering 
amplitudes. This diagram can be calculated by numerically integrating with the half off- 
shell LO scattering amplitude and has been done in Refs. 16|, |30|. At NNLO there are two 
diagrams one must calculate in Fig. HI The first diagram simply sandwiches two effective 
range corrections between two half off-shell LO scattering amplitudes and can again be 
solved straightforwardly numerically. However, the second diagram has an insertion of a full 
off-shell LO scattering amplitude and in principle cannot be solved without calculating the 
full off-shell LO scattering amplitude. In order to circumvent the need to calculate the full 
off-shell LO scattering amplitude the partial-resummation technique has been employed by 



Bedaque et al. 16j , which gives the correct diagrams up to NNLO but includes certain higher 
order diagrams and therefore is not strictly perturbative. A fully perturbative calculation has 
been carried out by Ji and Phillips j^j by effectively calculating the full off-shell scattering 
amplitude. However, their calculation was not for nd scattering and was only for S-wave. 
Now we will introduce a new approach to calculate the nd scattering amplitudes to NNLO 
strictly perturbatively that avoids calculating the full off-shell LO scattering amplitude. 
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FIG. 3: Diagrams at NLO 



TT 




FIG. 4: Diagrams at NNLO 
The full nd scattering amplitude is in general given by 

i, q (k,p) + t l hq (k,p) + t l 2tq (k,p) = B l (k,p) + B[(k,p) + B l 2 (k,p)+ (13) 
+ (K l (q,p, E) + K[(q,p, E) + K l 2 (q,p, E)) <g> (t l q (q, k) + t[ jq (q, k) + t^q, k)), 

where t l n (k,p) is the scattering amplitude, B l n (k,p) the inhomogeneous term, and K l n (q,p, E) 
the kernel of the integral equation. The subscripts 0, 1 and 2 refer to LO, NLO, and NNLO 
parts respectively. The symbol £g> represents an integration and is defined by 



A(q)®B(q) 



TT 



dqq 2 A{q)B{q). 



(14) 



The inhomogeneous part and kernel of the integral equation to each order in EFTy are 
defined by 



B l (k,p) 



y 2 t M N 

pk 



Qi 



p 2 + k 2 - M N E - it 
pk 



B l 1 (k,p) = 0,B l 2 (k,p) = (15) 



K{q,p,E) 



M N Vt n (n) 



DX 1 E 



_ \ 1 _ / q 2 + p 2 - M N E - it 



2M N J qp 



(16) 



47r " \ ~zivi N j qp \ pq 
At NLO in the partial-resummation technique one keeps the inhomogeneous part, kernel, 
and amplitude up to and including NLO. Denoting t l NLO = t l 0q + t[ q (Note for ease of 
discussion we have omitted the momentum dependence of all functions in the text) one 
obtains 



*WM = B l (k,p)+B[(k,p) + (K l (q,p,E) + K[(q,p,E))®t l NLO (k,q), (17) 



9 



for the NLO amplitude in the partial-resummation technique. This technique again is not 
strictly perturbative since the K[®t[ q term is NNLO. Thus in order to make the calculation 
strictly perturbative one simply throws out the term K[ <g> t[ q and obtains 

t[ q (k jP ) = B[(k,p) + K[(q,p,E) ®tl q (k,q) + K l (q,p,E) ®t[ jq (k,q), (18) 

for the NLO piece of the scattering amplitude. In this equation the term K[ <g> t Q „ is simply 
absorbed into the inhomogeneous part of an integral equation for t\ since the amplitude 
t l is already calculated at LO. The resulting kernel for the integral equation is K l Q , the 
same as the kernel for the LO amplitude. Now collecting all of the NNLO terms in Eq. ( TT5j) 
one obtains 



t[ ;q (k,p) = B l 2 (k,p) + K l 2 ® t^k, q) + K[(q,p, E) ® t l ^ q (k, q) + K l (q,p, E) <g> t l 2tq (k, q), (19) 

for the NNLO piece of the scattering amplitude. In this equation K 2 ® t l Q q and K\ £g> t[ q are 
absorbed into the inhomogeneous part of an integral equation for t\ , since t l is already 
calculated at LO and t\ is calculated in Eq ( !T8|) . Again the kernel for the integral equation 
of t\ q is given by K l . Thus at each order the kernel of the integral equation is always the 
same but the inhomogeneous part of the integral equation changes at each order and depends 
on lower order amplitudes. Therefore, in order to calculate the NLO and NNLO amplitudes 
numerically one can simply use the kernel calculated for the LO scattering amplitude as 
it is now the kernel for all orders. The Hetherington-Schick method can also be used to 
calculate the NLO and NNLO amplitudes. Since t l Q and K[ have no singularities in the 
contour used to perform the integration, the term K[ ® t introduces no new singularities 
inside the contour. Therefore, no new singularities are introduced in the contour for the t[ 
scattering amplitude, and the same analysis holds at NNLO. 

IV. Doublet Channel 

In the doublet channel two extra complications arise. One comes from the fact that we 
now need to include the singlet propagator, which gives two coupled integral equations. 
The second arises from three-body forces, which only occur in the doublet S-wave channel 
starting at LO because there is no centrifugal barrier and the Pauli exclusion principle 
does not exclude a bound state in this channel. Since the resulting kernel of the integral 

10 



equation in the doublet S-wave channel is non-compact, the equation does not possess a 
unique solution (in the limit A — > oo). At finite cutoff this results in the solution varying 
greatly with the choice of cutoff. This can be remedied by the insertion of a three-body 
force with an appropriate scale dependence on A. The basic form of the three-body force 



up to NNLO is given by [15|, |l7|, [31 1 



H(E, A) = °J } + ° A2 l ' + ° A2 1 ' + 2 A4 1 \ M N E + 7 2 ). (20) 

The term H is chosen at each order so that we get the correct scattering length of 
ai = .65 fm in the doublet S-wave channel. This parameter is split up into pieces depending 
on the order at which we are working, so that H doesn't have to be refit at each order. 
The value of H 2 is chosen such that we get the correct triton binding energy of = —8.48 
MeV. Unlike in the partial-resummation technique the binding energy must be calculated 

n 

perturbatively here, and we use the technique as outlined by Phillips and Ji [221 ] . Now with 
all these complications in mind we must solve a set of coupled integral equations represented 
in Fig. [5j Note that diagrams with three-body forces will only contribute to the doublet 
S-wave channel. 




FIG. 5: Integral Equations for doublet channel at LO. Note diagrams with three-body 
forces only occur in doublet S-wave channel 



In Fig. |5]the double dashed line represents the LO singlet propagator and the thick solid 
(blue) line represents a sum of two diagrams, one containing a deuteron propagator and the 
other having a singlet propagator in place of the solid (blue) line. 

11 



Having carried out all necessary projections in the doublet channel as in Refs. |16l . Il8| . 
it is convenient to represent the integral equations in cluster-configuration space (l8j. The 
master Eq. ( TT3"j) in the quartet channel now becomes (Note the subscript d on the amplitude 
t refers to the doublet channel.) 



tUM +t[ jd (k,p)+t l 2id (k,p) = B l Q (k,p) + B[(k,p) + B l 2 (k,p) (21) 
+ (K l (q,p, E) + K[(q,p, E) + K l 2 (q,p, E)) ® (t l 0jd (k, q) + t[ 4 (k, q) + t l 2>d (k, q)), 

where the bold script denotes that this is now a matrix equation in cluster-configuration 
space. The vector d (k, q) defined by 



t 



t 



n,Nt-*Nt 
I 

n,Nt^Ns 



(k,q) 



(22) 



contains an amplitude t l n Nt ^ Nt (k, q) (t l nNt _> Ns (k,q)) that corresponds to nd scattering (a 
neutron and deuteron scattering to a neutron and singlet dibaryon). The vector B n (k,p) is 
defined by 



B l (k,p) 



B[(k,p) 



Ay 2 M N 
■l2y t y s M_ 



n uQi 



r ^ ;{ ^-^- ) + n (E,A)5 m 



N 



pk 

p 2 +k 2 -M N E-ie 



pk 



+ n (E,A)6, 
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Ay 2 M N H l {E,A)8 l 
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-12y t y s M N n 1 (E,A)5 l0 
and the kernel matrix K. l n (k, q, E) defined as 



iy 2 M N H 2 {E,A)5 l0 
-12y t y s M N H 2 (E,A)5 l0 



(23) 



M NnM ,„ f q 2 +p 2 -M N E-ie 



rt n \EA)-Qi 

57r qp 



M 



qp 



-Vt ZytVs 
3y s yt -y 2 s 



+ ^my2^ U \EA)n n - j (E,A) 



2 



j=o \ yty s -y, 

where D^ n ^(i?, q) is a matrix of dibaryon propagators given by 



(24) 



A (n) (£,q) 



B n (E,q) 







Di n) (E,q) 



(25) 
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which multiplies the matrix occurring in the definitions of K n (k, q, E) via standard matrix 
multiplication to get the full form of K. l n (k, q, E). 

With the master Eq. (1211) in hand one can simply pick out all the strictly LO, NLO, and 
NNLO pieces to obtain an expression for the amplitude at each order, which are given by 



These equations are analogous to those of the quartet channel. The primary difference 
is that they are now defined in cluster-configuartion space. At each order the lower order 
amplitudes are put into the inhomogeneous part of the integral equation, and also the kernel 
at each order is always the same just as in the quartet case. 

Note that all of the amplitudes are unrenormalized in these equations. In order to find 
the renormalized amplitude we must multiply t l Nd ^ Nd (k,p) by the appropriate order of the 
deuteron wavefunction renormalization Zp. Lastly since the channel t l Nd ^ Nt (k,p) is of no 
direct physical interest to us we can renormalize it in any manner we wish. We choose 
to renormalize t l Nd-^m(k,p) such that after renormalizing t l Nd ^ Nd (k,p) by the appropriate 
order of Zd, all dependence on y s and y t are removed from our integral equations. 

V. SD mixing 

Finally at NNLO there is an SD-mixing term that will cause splittings in partial waves 
of different J values (J = L + S is orbital plus spin angular momentum) as well as mix- 
ing different partial waves with the same J values. The Lagrangian in the auxiliary field 
formalism for the NNLO SD-mixing term is given by 



t l 0jd {k, p) = B l (k, p) + K l Q (q, p, E) <g> t l 04 (k, q) 



(26) 



t[ jd (k, p) = Bi (k, p) + K[ (q, p, E) <g> t l 0>d (k, q) + K l (q, p, E) <g> t[/k, q) 



(27) 



H,d(k,p) = B l 2 (k,p) + K l 2 (q,p,E)®tl d (k,q) + 

+ K[(q,p, E) ® t\ d {k, q) + K l (q,p, E) ® t l 2 d (k, q) 



(28) 
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and in terms of nucleon fields the leading SD-mixing term in EFTy is given by 



w - ^(.Y'/^.Y) 



N 1 



(d - d)\d - d) J - -s i] (d - dY ) p 3 n 



+ h.c. 



(30) 



The parameter Cq°_ x is fit by producing the correct asymptotic D/S ratio in the deuteron 
wavefunction Q which yields 



r<SD 
°o,-i 



-Vsd- 



(31) 



l M Nl ^-ltY 

where 77^ = .02543±. 00007 [32] is the asymptotic D/S mixing ratio of the deuteron wavefunc- 
tion. Integrating out the auxiliary fields and performing a field redefinition of the nucleon 
field one can relate i/sd to 1 via 



1 A (3Sl) 

VSD _ 1 ^(-1) riSD 
,2 °0,-l- 



(32) 



Vt 4 yf 

Then using Eq (IHTj) one can obtain a value for the parameter ysolvt which will ultimately 
appear in all of our SD-mixing amplitudes as we will show below. 

The resulting amplitude from the SD-mixing term is given by the sum of diagrams in Fig. 
EJ Note that these diagrams are not projected out in either angular momentum or spin, and 
upon projection in certain channels specific diagrams will give zero contribution. 





FIG. 6: SD-mixing diagrams at NNLO (Note diagrams where lower vertices contain 
SD-mixing terms are not shown) 



The sum of all diagrams in Fig. [6] in cluster-configuration space gives 
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where the Greek letters represent nucleoli spinor indices, Latin letters a,b,c, and d the 
nucleon isospinor indices, and finally Latin letters x,w,y, and z the deuteron or dibaryon 
polarization depending on the respective matrix element in cluster-configuration space. As 
before, the vector k (p) represents the the momentum of the incoming (outgoing) particles 
in the cm. frame. The vector v p projects out the nd amplitude in cluster-configuration 
space and is defined by jisj 




Vp= l v oJ' (34) 

and the amplitude y(it xw )^ a (k, q)J is a vector defined by 



(A+XW \P b 
1 LL Nd-^Nd)aa 




where t Nd ^ Nd is the amplitude for nd scattering and t^d^m is the amplitude for nd going 
to a nucleon and a singlet combination of remaining nucleons. (Note that we have not yet 
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projected out quartet or doublet channels.) The matrix D(°>(E, q) was given earlier in Eq. 
(1251) . and (iK )^(q, £, go) is a matrix defined by 



X 



' aa 



(36) 



where the functions (JCgp xw ) aa (q, £), which contain all of the angular dependence are de- 
fined as 



>a u a 



(2p + k) w (2p + k)» - i^(2p + k) 2 



(2k + pf (2k + p)« - -5 yx (2k + p) 2 
(2k + p)* (2k + p) y - ^6 yx (2k + p) 2 
(2p + k)-(2p + kf -I^(2p + k) 2 



(37a) 

(37b) 
(37c) 
(37d) 



(Note that the capital letters A, B, and C are used for the singlet auxiliary field polarization 
and the lowercase letters w,x, and y are used for the deuteron auxiliary field polarization.) 
Integrating over the energy and picking up the poles from the nucleon propagators in our 
diagrams, Eq. (|33|) becomes. 



«> ,C ^ _ M N T (K xw 

2 p \ /aa 



M 



+ 



k,p)v p (38) 
2 / ^(^m P)D (0) (E - q) ((f)- (k, q 

M " ' ^ ((t-)J(q,p)) T D(°)^- ^,q)(/CnS(k,q)v P 

d 3 / / „, _ xT 



jV 



2 



(2tt) 3 
d 3 g 



/(S/(S(rt^')) TD(o, (-4-) 

(/C^)^ c (q,£)D(°) (^-^-^) ((^XaAq 
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where 



(K xw )ia(<iJ) = ^ x (39) 

q 2 + q_-£ + £ -M N E-ie 

x ( m l D XW )l(^) {K>f D ™fl (q^) \ 

X V {K%r)£ (q^) (4?) / ' 

Now since the SD-mixing term mixes spin and orbital angular momentum we will project 
out in total angular momentum J, and the amplitude for the SD-mixing term can be written 
in a partial wave basis as 

oo M=J J+S J+S' 

tsD (k,p) = j2J2 E E E 4 ^^^L5(^p)^(p)(^/^ s (k))* (40) 

J=0 M=-J L=\J-S\ L'=\J-S'\ S,S' 

where the spin angle functions are 

0#s(k) = £ C™' M lT L (k)xr- (41) 

Here, Xs s * s the spinor part of the spin-angle functions, C™s'™ s,M the appropriate Clebsch- 
Gordan coefficient, and Y™ L (k) the appropriate spherical harmonic. Since the spin-angle 
functions are orthogonal, we can project out the amplitudes in our angular momentum basis 

tsDi&jjfap) = ^JdQ k j dQ p {&j M L , s ,(p)yt SD (k, p)^| 5 (k) (42) 

Finally projecting out the isospin and carrying out the spin angle projection on Eq. fl3"5|) 
we obtain 



WfiWM = ^ T p m,P)-hs',LS^+ (43) 
- 1^ jf dqq2 ^,L'S^P)) T ^ ] ( E ~ <^,^ 
x /C(g,£)i*l, LS D(°) [ £ - * 2) (t™ , 
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where t^| f i5 (/c, q) are the LO scattering amplitudes projected out in a partial wave basis 
and have already been calculated in the preceding sections. The matrix K,(k,p) J L i SI LS is 
the kernel function of Eq. (l3T?j) fully projected out in isospin and total angular momentum. 
The kernel matrix (/C a:w )^(k, p) can be fully projected out in isospin and total angular 
momentum by use of Racah algebra obtaining general expressions for any values of angular 
momentum in terms of 3n-j symbols and has been done in detail in Ref. 



. The matrix 



element [7C(&,.p)ls,l'S'] 2 2 = ® ^ e subscript outside the square brackets represents the 
specific matrix element) since there is no two-body SD-mixing term coupling to the ^So 
dibaryon. The matrix elements [/C(/c,p)£ S)i / S /] and [/C(fc,p)^ Si , 5 ,] 21 are related by 

M*, P)isM 12 = [*(p. k)l s ,, LS ] 21 , (44) 

because of time reversal invariance. Finally \fc-(k,p) J LSUS ^ and \K-{k,p) J LStL i S i\ are given 
by 



[K(k,P)t> S >, L s]u = 4 7 rV / M^ y (<V,y 2 + 2S s ,, 3/2 ) C^i-l) 28 '" 

2 1 1 1 I S' 2 S I 1 

x< }{ \]Tr(k 2 Qu(a)+4p 2 Q L (a)) 

1/2 5 S' \ L J V \ k P 



-S+L-.J 



X 



+ 8tt y/B&L (6 s ,,i /a + 28 s , j3/2 ) C L°i>> C L»iL' 



x 



X 



1/2 1 S 

1 L" L 
S' V J 



> Ql"{cl) 



x VL 



1/2 1 S' 



L 1 L" 



V J L" 1/2 J S 



Ql"(o) 



16tt 1 



(<V,y 2 + 2S 5 ,, 3/2 ) (-l) 1/2 - S '5 L , L ,5 s ,s> J2 C LZ>> C lL»,LVL"Q L >>(a) 



3 VI V S ' 1/2 
+ (5 ^ — y S'){L < — y L'){k <— ► p) 



(45) 



and 



18 



[K(k,P) J L'S',Ls} 21 = SnV5VSL6 S n /2 C°^%(-l) 



,1+S+L-J 



(46) 



2 1 1 5" 2 S 1 , 
1/2 5 5' J \ L J L' \ k P 

* ' 1/2 1 5 

1 L" L 

v 

tor>/Sy/LS6a,y* ^(-l^+^/^W r ,o.n,, 



> Q L »(a) 



^ °L,1,L"°L",1,L' X 



L' 1 L" 



L 1 L" 



Ql"{o) 



167T , lN r// r 

+ vf M) vi 



1/2 J S" I I 1/2 J S 

F7 V ^^L,i',L''^L",1,L'^'V2^5'5^L'L<5l 



where Ql(gi) are defined by Eq. ([12]) and a = (A; 2 +p 2 — M N E)/pk. 

In order to evaluate Eq. (1431) we multiply both sides by Zlo, and use the renormalized 
scattering amplitudes in the integrals. Doing this this will remove all factors of y s and all 
terms will be multiplied by factors ysD/yt, which has already been fit to physical data. 
Then numerically integrating the renormalized scattering amplitudes with their matching 
kernel we obtain the renormalized SD-mixing amplitude. Note that the renormalization of 
the t l n Nt ^ Ns (k,p) was arbitrary and chosen to remove any dependence on y t and y s in the 
doublet integral equations. As long as we are consistent with our normalization choice all 
factors of y s should also cancel out in the SD-mixing terms leaving only yso/yt 

VI. Phase Shift Analysis 

After calculating the quartet, doublet, and SD-mixing amplitudes one can calculate the 
scattering amplitude T[, s , LS up to NNLO, where L (S) is the initial orbital (spin) angular 
momentum, L' (S') the final orbital (spin) angular momentum, and J the total angular 
momentum. At LO the scattering amplitude is given by 

t usIls(p) =hL'Sss'Ss'y2Z LO to, Nt ^ Nt (p) + S LLI 8 SS '5 S '3/ 2 Z LO to iq (p), (47) 
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where the quartet and doublet scattering amplitudes are multiplied by the LO deuteron 
wavefunction renormalization. (Note the superscript (n) on T^jg! LS refers to the order of 
the amplitude.) The NLO scattering amplitude is given by 



t l'S',ls(p) =hL'5ss'Ssn/ 2 (Z NLO tl Nt ^ Nt (p) + Z LO t^ m ^ m {p))+ (48) 
+ S LU 5 ss/ 5 S '3/ 2 {Z NL oto !q (p) + Z LO t\ q {p)). 

For both LO and NLO the scattering amplitude has no J dependence and there is also no 
mixing between different spin and partial waves. Finally the NNLO scattering amplitude is 
given by 



t l'S',ls(p) =hL>8ss>8s'i/2(Z NNLO t^ m ^ Nt (p) + Z NLO t^ Nt ^ m {p) + Z LO t^ Nt _> m {p))+ (49) 
+ hL'Sss'S S '3/2(Z NNLO t^ q (p) + Z NLO t{ q {p) + Z LO t^ q {p)) + t SD i, s , >LS (p), 

and at this order the SD-mixing amplitudes are necessary. Note the SD-mixing amplitudes 
in Eq. (|4"9"1) are already renormalized by the LO deuteron wavefunction renormalization 
and thus do not have a factor of Z^q. The SD-mixing amplitude at this order will cause 
splittings for different J values in the quartet channel. However, the SD-mixing amplitude is 
zero for the doublet channel and no splitting occurs. Also the SD-mixing amplitude allows 
all possible (up to conservation of parity and total angular momentum) mixings of different 
initial and final partial waves and spin. 

Since the total angular momentum is a good quantum number the S-matrix can be 
decomposed into a direct sum of irreducible representations of J. For each J value, states 
of different parity can be further factorized since parity is conserved. Starting at J = | the 
S-matrix is decomposed into a 2x2 matrix for positive and negative parity as in Eq. ( 150]) . 

(Q2 Q2 \ / Q2 Q2 

ft p >S*"= V* I (50) 

I I 

The 2x2 S-matrix in the Blatt and Biedenharn parametrization [34j is given by 

S J -= (u^)V i<5 '%^ (51) 
where the 5 Jn and u j7r (note tt represents parity) are matrices defined as 
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•< I ! 5 .j I » I 'J ^ I. (52) 

i , / cosr?^ + sinr?^ + \ i / cose^~ sine^" 
\ — sinij2 + cos 7/ 2 " + " / \ —sine 2 cos 6 2 

In the limit where the mixing parameters and are zero the matrix u j7r becomes 
the identity matrix and the S-matrix becomes diagonal. Therefore, the 5 J L n s eigenphase 
parameters become the standard phase-shifts in the limit of zero mixing angles. To obtain 
fits for the phase-shift parameters one must relate the S-matrix to the T-matrix and this 
can be done by projecting out the operator equation S = 1 + iT in total angular momentum 
yielding 

Sl's>,ls = Sll'Sss' + i^-^-Ti, s , LS (p). (54) 

Now one expands the eigenphase parameters 8j]g and the mixing parameters perturbatively 
up to NNLO giving the S-matrix in a perturbative expansion. Then using Eq. (1541 one 
expands the S-matrix perturbatively by expanding the T-matrix perturbatively and making 
5ll>5ss' LO. Finally matching these two perturbative expansions onto each other one obtains 



^» = ilog(l + ^rS s (p)) (55a) 



1 ■2M M p T J{l) / \ 
f J(l) _ 1 1 3tt 1 ls,ls\P) 



3tt ^W-S 

for the perturbative expansion of the eigenphase parameters and 

T V2(2) T V2(2) 
Z 2' Z 2 U 2' U 2 i 2' i 2 i 2' i 2 

for the mixing parameters, in terms of the T-matrix. Note that the mixing parameters first 
occur at order NNLO. 
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For J > | the S-matrix is given by a 3x3 matrix and we choose to parametrize the matrix 



using the conventions of Seyler [35], in which case the S-matrix takes the form 



D Jin3 3 r-p33 
J i 22' +22 



s 



■it 



2 2 2 2 



11 S T . 



■IT 



2 2 ,J=I= 2 2 



S 



\ J=c 2 2 



11 '-'Till 



5 



1 1 S J T ,i 



(57) 



To parametrize the matrix we again use Eq. (1511) . However, the matrix S Jn is now a diagonal 
3x3 matrix defined by 



/ 5i 



;Jtt 



V 



o 
o 



o 

J=c 2 2 

5 






1 



\ 



J 



(58) 



and we now have three mixing parameters e Jn ,( Jn , and rj Jn . The matrix u Jn is now defined 
by three successive Blatt and Biedenharn "rotations" 



cosC j7r sin( j7T \ 

V Jn = I cose j7r sine^ , W J?T = I 10 , X. Jn = I -sin^ cost?- 711 I (59) 

-sint j7T cose j7T J \ -sin( j7r cos( j7T J \ 1 

u J " = v J "w j7r x j7r (60) 

Finally, expanding the 3x3 S-matrix perturbatively in terms of the phase-shift parameters 
and the T-matrix we can again obtain an expression for the perturbative expansion of the 
phase-shift parameters in terms of a perturbative expansion of the T-matrix. The equations 
for the eigenphase parameters are the same as before. However, the mixing parameters are 
given in terms of the T-matrix as follows 



J +22 



J=c 2 2 



T J(2) 



' J=c 2 2 



T J(0) 

1 Jzp33 / T 33 
J T 2 2' > 2 2 



T J(0) 

' — T I 1 ' 



2 2 ' 2 2 



JX 22 ' 2 2 

rp/(0) _ j,J(0) 

"22' 22 22 



7^(2) 

J JxM t ± 13 

J '22> 22 



T J(0) 



3 3 



J=P^3 
2 2' J_I "2 2 



2 2'" y '-'- 2 2 



(61) 



JI 22 
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VII. Results 

The only channel for which three-body forces are necessary up to and including NNLO 
is the doublet S-wave. After fixing the three-body forces as described above the resulting 
doublet S-wave phase-shift is given in Fig The data below the deuteron breakup thresh- 



old DBT is from the AV18+UIX potential 



with the Pair Correlated Hyperspherical 



Harmonic method (crosses ), a nd above and below breakup from the Bonn-B potential with 
Faddeev equations (stars) 37|]. Agreement between the potential model calculations (PMC) 
and the calculation here is within the expected error of EFTy in the Z-parametrization at 
NNLO, {{Z t - l)/2) 3 « 3%. 
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FIG. 7: (Color online) quartet (doublet) S-wave phase shift above (below), real part 
(imaginary part) right (left) in degrees as function of lab energy, quartet S-wave given for 
A = 1600 MeV and doublet S-wave cutoff varied from A = 200 — 1600 MeV. The crosses 
below breakup are AV18+UIX data 36], the stars Faddeev equations with Bonn-B [3 
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At each order we vary the cutoff A between 200 and 1600 MeV. The resulting range of 
predictions is shown as a solid band in the figure. One should also note that each order 
the cutoff variation diminishes. Also shown in Fig [7] is the quartet S-wave phase shift 
which agrees well with the same PMC. The cutoff variation in the quartet S-wave channel is 
negligible as no three-body forces are necessary, thus the phase shift is plotted for a cutoff 
of A = 1600 MeV. Note at roughly A = 900 MeV all results begin to converge even in the 



doublet S-wave and higher partial waves, which is in agreement with previous findings [18]. 
The scattering length, 4 a, in the quartet S-wave channel at NLO and NNLO is found to 
be 6.74 fm and 6.19 ± .030 fm respectively. The NNLO scattering length compared to the 
experimental value of 6.35 ± .02 fm 38) is within the expected 3% error. In the partial- 



resummation technique the NNLO scattering length of 6.354 ± .002 fm [18J] is much closer 
to the experimental value and is not surprising since the partial-resummation technique 
contains certain higher order contributions. 

For higher partial waves in the quartet channel there will be splittings for different J 
values due to the inclusion of the SD-mixing term at NNLO. In Fig. [8] the eigenphases for 
P through G-waves are shown for each possible J value and are compared to data from 
AV18+UIX below the DBT 36(. Agreement between the PMC is apparent for partial waves 
D through G both qualitatively and quantitatively. The AV18+UIX data agree within the 
predicted 3% error and the direction of the splittings is also consistent with our calculation. 
However, for the P-wave there are notable discrepancies, namely the AV18+UIX data does 
not agree within the predicted 3% for larger energies below the DBT or put simply our 
splitting appears too large. Despite the large splittings there is still qualitative agreement in 
the direction of the splittings for different J values. For Fig [9] the eigenphases for P thr oug h 



G-waves are shown for energies above the DBT with data from the Bonn-B potential |37 |. 

The qualitative agreement between the PMC and our calculation at these higher energies 
is worse since the splittings in the PMC cross each other and we find no such effect. At the 
higher energies the quantitative agreement with the Bonn-B data gets poorer (not within 
3%) as one approaches the breakdown scale of EFTy. One curious feature of the Bonn-B data 
is the J = | imaginary part of the F-wave eigenphase is much larger than all of the other J 
values for the imaginary part of the F-wave eigenphases. This is in clear contradiction with 
our results that show the splittings to all be of roughly the same size at this order in EFTy. 
Note that the issues in the P-wave noted below the DBT are only further exacerbated above 
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the DBT. 

Higher partial waves for the doublet channel are shown in Fig. (TDJ At NNLO the SD- 
mixing terms do not cause splittings in the doublet channel thus the data from AV18+UIX 

3?| above and below DBT are averaged for all J values to 



36[ below DBT and Bonn-B 



compare to our results. For the D through G partial waves good agreement is seen both 
for the real and imaginary parts within the expected 3% accuracy. However for the doublet 
3 -wave there are apparent discrepancies above the DBT in agreement with previous results 
1^ | . This may be due to the fact that the rather strong two-body P-wave contributions 
occur at N 3 LO and their inclusion may resolve this discrepancy. In addition this could 
explain the discrepancies observed in the splittings of the quartet P-waves. Thus further 
study at higher orders is warranted. 

We now consider the mixing angles, starting with the 7] J+ shown in Fig. [TTJ, which 
mixes states with the same partial waves but different total spin. The cutoff variation (from 
A = 200 MeV to A = 1600 MeV) for the r]i + mixing angle is shown explicitly since it is 
the only mixing angle that depends on the doublet S-wave. With the exception of the 77 1 + 
mixing angle good qualitative agreement with PMC can be seen for the real part of the t] J+ 
mixing angle and at low energies the agreement is within 3%. However, for the imaginary 
part of r] J+ at higher energies agreement with PMC quickly breaks down. In Fig. [Tj2] the 
mixing angles rj J ~ are shown and again good qualitative agreement is seen with the PMC, 
and at low energies the agreement with PMC is within 3%. However, at higher energies the 
discrepancies seem to be more pronounced than in the r] J+ mixing angles. This is likely due 
to the fact that these mixing angles depend on the P-waves and as already shown there are 
issues with P-waves that will hopefully be resolved with the inclusion of two-body P-wave 
terms at N 3 LO in EFTy. 

The mixing angles ( J+ shown in Fig. [13] mix states of different angular momentum but 
the same total spin. Good qualitative agreement for the real part of the mixing parameters 
C J+ is seen with the PMC and at low energies the agreement is within 3%. However, the 
imaginary part of £ J+ is roughly an order of magnitude smaller than the Bonn-B data at 
higher energies. In Fig the mixing angles ( J ~ are shown. For the real part of the 
agreement with PMC at low energies is apparent and at higher energies this breaks down, 

5 

in addition the imaginary part of has an opposite sign compared to the PMC. For the 

3 

real part of the qualitative dip seen in the potential model is not followed. However at 
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FIG. 8: (Color online) J splittings in higher partial waves below DBT. The dotted line is 
I — I, dashed line I — |, dot-dashed line I + 1, and solid line I + |. The points from 

36] are cross is Z — §, star I — |, open square J + |, and open triangle / + §. 



2 ' 

AV18+UIX 
A = 1600 MeV 



low energies there is still rough agreement. The imaginary part of £2" like £2 + seems to 
be an order of magnitude smaller than the PMC. Analogously to the 7] J+ and r] J ~ mixing 
angles it appears that the real part of ( J+ matches well to potential model calculation while 
( J ~ is qualitatively worse. This again is likely due the fact that ( J ~ depends on P-waves. 

Finally we examine the e J+ mixing angles in Fig. [15] as well the mixing angles e J ~ in Fig. 
[16] which mix different partial waves but the same total spin. For both sets of mixing angles 
agreement with the PMC is poor and the qualitative behavior of the PMC is generally not 
followed. Even at the lowest energies there is poor agreement with the PMC and the PMC 
are either over or under predicted. There is much need for improvement in these mixing 
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parameters and thus a higher order calculation of these mixing parameters is necessary. An 
accurate determination of all mixing angles and eigenphase shifts will allow the prediction 
of polarization observables in nd scattering and in particular the A y asymmetry. 

VIII. Conclusion 

The technique outlined here has general applicability in including perturbative corrections 
to integral equations. It has all of the numerical savings of the partial-resummation technique 
but with the added benefit that it is strictly perturbative. The partial waves obtained using 
this technique are very similar to those produced in the partial-resummation technique and 
similar issues are found in the doublet P-wave. As was shown, the quartet scattering length 
in the partial-resummation technique is closer to the experimental value at NNLO than ours. 
However, this is not surprising since the partial-resummation technique does contain certain 
higher order terms. More importantly the imaginary part of the quartet S-wave phase shift 
has the correct sign at NLO and NNLO in this technique, unlike the partial-resummation 
technique. In order to obtain the correct sign in the partial-resummation technique one 
must resum the effective range to all orders in the deuteron propagator and this introduces 
spurious poles which complicates the numerics. Note in this technique the SD-mixing terms 
could have been included in the integral equation for the NNLO amplitude and will be 
in future when higher orders are calculated. Also this technique allows one to calculate 
diagrams that contain external currents with full off-shell scattering amplitudes, thus in 
principle enabling one to calculate the three-body process 3 ife + 7 — >■ p + d , 3 H + ^y — >■ p + d, 
or even Compton scattering off of 3 H or 3 He. In short this technique makes the perturbative 
calculation in three-body systems numerically simpler. 

The calculation of nd scattering up to NNLO including the SD-mixing terms yields good 
agreement with PMC for the eigenphase shifts. However, as noted in the quartet P-wave 
the splittings are larger than expected and in the doublet P-wave where no splittings occur 
at this order in EFTy the potential models diverge from our results above the DBT. These 
discrepancies are likely due to the absence of the rather strong two-body P-wave terms that 
occur at N 3 LO in EFTy. Thus a higher order calculation to resolve this issue is in order. 
For the real part of the mixing parameters rj j7T and ( j7T , good agreement with the PMC 
at low energy has also been shown. As for the imaginary part of r/ j7T and ( j7r there are 
noticeable discrepancies with the PMC, in particular the imaginary part of ( j7r seems to be 
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an order of magnitude smaller than the potential model prediction. Also for negative parity 
r] J ~ and there is more noticeable discrepancy at higher energies with PMC and this is 
again likely related to the fact that these mixing angles depend on P-waves, thus further 
warranting the inclusion of higher order two-body P-wave corrections in EFT^. Finally the 
mixing parameters e Jn seem to match very poorly to the potential model calculation and 
could potentially benefit the most from a higher order calculation. 

In order to calculate polarization observables it is important to accurately determine the 
mixing parameters since in the absence of these there are no polarization observables. The 
technique outlined here will ease the calculation of higher order contributions necessary for 
an accurate determination of the mixing parameters. Of particular interest is to carry out 
a calculation for the mixing parameters in pd scattering as there are far more experimen- 
tal measurements to compare with. However, this will be complicated by the inclusion of 
Coulomb interactions. The technique of Hammer and Konig |20j can be combined straight- 
forwardly with the technique in this paper to calculate pd scattering strictly perturbatively 
in both Coulomb interactions and EFTy. However, at low energies it may be necessary to 
include the Coulomb interactions nonperturbatively and such a calculation is complicated 
by the numerical singularities introduced by the full off-shell Coulomb propagator, and new 
numerical techniques will need to be developed in order to deal with it. 
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FIG. 9: (Color online) J splittings in higher partial waves below DBT. The dotted line is 
I — |, dashed line I — |, dot-dashed line Z + |, and solid line Z + |. The points from Bonn-B 



37j | are cross is I — |, star Z — |, open square Z + 1, and open triangle I + |. A = 1600 MeV 
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FIG. 10: (Color online)Higher partial wave doublet phase shifts as function of lab energy 
real part (imaginary part) right (left) cutoff of A = 1600 MeV. Data shown below DBT is 
AV18+UIX (crosses) [36], above and below DBT is Bonn-B (stars) 37]. Note average 
((I + J) + (I - J))/2 of data is taken 
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FIG. 11: All rj J+ for A = 1600 MeV right (left) real part (imaginary part). 772+ has cutoff 
variation from A = 200 - 1600 MeV. Data below DBT (crosses) is AV18+UIX above 
and below DBT (stars) is Bonn-B 37] 
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FIG. 12: All 7] J ~ for A = 1600 MeV right (left) real part (imaginary part). Data below 
DBT (crosses) is AV18+UIX Q, above and below DBT (stars) is Bonn-B Q 
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FIG. 13: All C J+ for A = 1600 
DBT (crosses) is AV18+UIX 



MeV right (left) real part (imaginary part). Data below 
36| . above and below DBT (stars) is Bonn-B [37| 
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FIG. 15: All e J+ for A = 1600 
DBT (crosses) is AV18+UIX 



MeV right (left) real part (imaginary part). Data below 
above and below DBT (stars) is Bonn-B [37 1 
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FIG. 16: All e J - for A = 1600 
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